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Abstract 

Wc derive global analytic representations of fundamental solutions 
for a class of linear parabolic systems with full coupling of first or- 
der derivative terms where coefficients may depend on space and time. 
Pointwise convergence of the global analytic expansion is proved. This 
leads to analytic representations of solutions of initial-boundary prob- 
lems of first and second type in terms of convolution integrals or con- 
volution integrals and linear integral equations. The results have both 
analytical and numerical impact. Analytically, our representations of 
fundamental solutions of coupled parabolic systems may be used to 
define generalized stochastic processes. Moreover, some classical an- 
alytical results based on a priori estimates of elliptic equations are 
a simple corollary of our main result. Numerically, accurate, stable 
and efficient schemes for computation and error estimates in strong 
norms can be obtained for a considerable class of Cauchy- and initial- 
boundary problems of parabolic type. Important instances of applica- 
tion are representations of solutions of multidimensional Burgers equa- 
tions with forcing and potential initial conditions and Pauli equation 
describing the non-relativistic limit of Dirac theory for electrons in a 
magnetic field. 

Warning: The current analysis holds only in special cases, essen- 
tially in the case of scalar equations. A more involved analytic expan- 
sion is necessary (and possible) for systems and will be communicated 
soon. 
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1 Introduction 

We consider linear equations of the form 
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^ = V 2 u + B.Vu (1) 

on a domain D = (0, T] x Q with $7 C R n a bounded domain, and where 

u = (-ui, • • • ,u n ) T (2) 

is a vector-valued function and B = (B , • • • , B n ) is an n-tuple of matrix- 
valued functions B l = (6*- fe ) where each entry b % - k possibly depends on space 
and time. More precisely, we understand B • Vu as a vector the ith compo- 
nent of which is given by 

n ft 

I™ (3) 

j,k = l 

such that general linear coupling of first order terms can be expressed. This 
means that in coordinates equation (pQ) is given by 

at dx* + jk dx k [ > 

3=1 ^ j,k=l 

for 1 < i < n. We are looking for an analytic representation of the solution 
(t, x, s, y) — > p(t, x; s, y) for ([1]), (J4j) with Dirac distributions 8 y {x) = 5(x — y) 
as initial conditions, i.e. for a representation of the fundamental solution. In 
the time- homogenous case dependence of time is only dependence of t — s, 
so that in this case we also write the fundamental solution in the form 
(t, x, y) — > p(t, x; y). For our representations of the fundamental solution for 
equations with time-dependent coefficients we also fix the parameter s = 
and write the fundamental solution in the form (t, x, 0, y) — > p(t, x; 0, y) for 
simplicity of notation. 

Remark 1.1. In the following we denote the fundamental solution of a system 
by bold face letters and use usual type of letters for the fundamental solution 
of a scalar equation. 

Remark 1.2. In general for parabolic systems of order 2p of form 

8 N 

^=EE m^)9>3 (5) 

j=l \a\<2p 

(with some natural number iV) the fundamental solution (or fundamental 
matrix) (t, x, s, y) — > p(t, x, s,y) is a N x A^-matrix of functions on (0, T] x 
O x (0,T] x CI which satisfies ([5]), and such that 

lim / f(y)p(t,x,s,y)dy = f(x) (6) 



2 



for all continuous functions y — ► f(y) in f2. The latter condition is equivalent 
to the rule that p(0, x; 0, y) = 5 y (x) = 5{x — y). Here in the general case with 
coupling of the higher order derivatives 2p, p > 1, a vectorial representation 
of the fundamental solution is not known. It is one of the observations of 
this paper that a vectorial representation is possible in the case p = 1, if the 
only coupling occurring is that via first order terms. Note that N ^ n is 
possible. Our restriction to the case iV = n is not essential but only related 
to an economy in the use of symbols. 

Remark 1.3. Probabilistic representations of the solution of the Cauchy 
problem for linear parabolic systems have been obtained first by Stroock 
(as far as I know, cf. [13] )• Note that in [J3] a representation is obtained via 
a fundamental matrix representation different to our representation which 
is vectorial. Moreover, our approach may be used in order to extend the 
Feynman-Kac formalism slightly beyond the class of processes considered 
in [13], because we may have n diffusion matrices a l - k , 1 < i < n in our 
generalization mentioned in the last section below. 

It turns out that results in the case of time-homogenous coefficients can 
be extended to the case of time-inhomogeneous coefficients but it is worth 
to consider the time-homogeneous case separately, because less assumptions 
have to be made. For this reason we shall state our main theorem in the 
time-homogenous case separately, i.e. where coefficient functions are of the 
form x — ► bj k (x). We shall assume that the functions x — ► b^ k {x) and their 
derivatives are uniformly bounded by powers of a generic constant c such 
that 

\d«b) k \ < C H (7) 

for all multiindices a = (ai,--- ,a n ). Here = Q x ^i^ dx a„ denotes the 

partial derivative operator of order a with respect to x = (x±, ■ ■ ■ ,x n ). If 
the coefficients are time-dependent functions of form (t, x) — > bj k (t,x), then 
we shall assume in addition that 

\dTb) k \ < Cm\ for all integers m > 0, (8) 

all 1 < i,j,k < n, d™ = J^r is the derivative of order m with respect to 
time. 

Note that assumption [7] holds for arbitrary finite Fourier series. Clearly 
it holds also on a bounded domain for polynomial coefficients, and on such 
domains multivariate polynomials can approximate all continuous functions 
in the supremum norm. This means that we are flexible enough for numerical 
applications. Indeed the treatment for higher dimensional scalar parabolic 
problems based on analytic expansions of the type considered here (simpli- 
fied to the scalar case) showed accurate and fast computations as well as 
error estimates in strong norms (cf. [3]). The main reasons, however, to 
introduce assumption [7] are analytical. First assumption [7J implies that b % - k 
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are globally analytic, i.e. for all y £ M. n b l - k equals its Taylor expansion, i.e. 
we have 

b) k {x) = £ %^A*«, (9) 
J L — ' a! 

|q|>0 



where a is a multiindex and Ax = (x — y), and bj ka (y) = d^bj k (y). This 
makes it possible to write down explicit solutions of parabolic systems of 
typed] in terms of power series of type El Second, the proof indicates that 
the assumptions made here cannot be weakened in general. 

Essential parts of our considerations can be generalized to parabolic sys- 
tems with space-time dependent second-order terms, i.e. equations of the 
form 

^i = f „j t ^L + f ij^, (io) 
dt ^ 3k dxjdx k ^ 3 dx k 

j,k=l J j,k=l 

where the scalar functions a l - k may depend on space and time. It turns 
out that the convergent analytical solutions (in case ([!])) are building blocks 
for the representations of solutions of parabolic systems with potential and 
source terms. As examples, let us consider a Cauchy problem and a standard 
initial-boundary problem which occur in the vector-valued as well as in the 
scalar case. The finite horizon Cauchy problem for parabolic systems of type 
© is 



& " E U & " E jk b} k ^ = fi m R» x (0, T] 



(11) 



_ Ui(0,x) = (j>i(x) 
where T > and for 1 < i < n. 



on 



Remark 1.4. The proof of the pointwise valid representation of the funda- 
mental solution given is valid for bounded domains Q, and cannot be directly 
generalized to unbounded domains. This is no essential restriction for nu- 
merical treatment, however. Analytically, a generalization is possible, if one 
considers a slightly different representation (cf. section 4.2.). However, the 
recursive relations of the expansion coefficients are more complicated and 
the convergence proof is more involved. Therefore we restrict ourselves to 
the case of bounded domains 0, in this paper. 

Another example is the initial-boundary problem of second type. We 
consider it in the scalar case here. Consider a domain C W 1 and denote the 
three constituents of boundary of the cylinder by := {(t, x)\t = & x € Q} 
Ox(0,T) by T = {{t,x)\t = T & x £ Q} and 5 = {(t,x)\t £ (0,T) & x £ dfl}, 
where 90 denotes the boundary of 0. The initial-boundary problem is of 
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the form 



( du 




n d 2 u 



Eu du 



f 



in n x (0, T] 



«(o, .) = 0(O 



on 0, 



(12) 



v dt 



on B 



where a, <ft, and ^ may depend on space and time. With an explicit repre- 
sentation of the fundamental solution we can represent the solution of ([lip 
in terms of convolutions of the initial data and the source data with the 
fundamental solution, and the solution of (|12j) in terms of convolutions of 
initial data, source data, and a function which is solution of a linear inte- 
gral equation. It is clear that such representations lead to accurate schemes 
which have obvious advantages compared to finite difference schemes and 
other standard schemes. 

Remark 1.5. In ([!]) we may add potential terms of form CiU with a coefficient 
functions c; which may depend on space and time. Theorem 1 below can be 
trivially extended to this case. Hence, in equations (jlip and ()12p we may 
also add potential terms of form c^u and representations of solutions in terms 
of convolutions and linear integral equations (in case of the initial-boundary 
problem (fT2j) can be obtained. 

Remark 1.6. More general cylinder domains D = Uo<t<T^t with varying Qt 
may be considered, of course. 

This is the first paper on globally pointwise valid analytic expansions 
of parabolic systems. In the case of scalar equations there are some in- 
vestigations and applications to problems of computation recently ([3] and 
references). Our result has direct applications to case of the scalar equa- 
tions, of course. Further comments on the relation to results in the scalar 
case can be found in Section 7. 

The outline of this paper is as follows. In the next section we state the 
main results concerning the representation of the fundamental solution. In 
Section 3 we formally compute the analytic expansion of the solution and in 
Section 4 we prove the pointwise convergence of the analytic representation 
in the time-homogenous case for a certain limited time horizon < t < Tq. 
In Section 5 we extend the results of the preceding Sections to the case where 
the coefficients may depend on space and time and we show the global con- 
vergence for any time horizon < T < oo. In Section 6 we consider the 
implications for representations of solutions Cauchy problems and second 
initial- value boundary problem and briefly discuss the advantages for build- 
ing efficient numerical schemes. In Section 7 we state some generalizations 
with general (but uncoupled) diffusion coefficients and show that a result by 
Varadhan is a rather immediate consequence of our main theorem. We also 
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discuss possible other applications (for example the definition of generalized 
processes) and give some further comments and an outlook. 



2 Main results on linear parabolic systems 

Since the second order derivative terms in Q are uncoupled, we may expect 
that a vectorial representation of the fundamental solution p = (pi, ■ ■ ■ ,p n ) 
(instead of an n x n fundamental matrix) is possible. The natural candidate 
for such a representation (in the time-homogeneous case) is 

1 / V n At 2 00 \ 

Pj (t,x,y) = ^rexp L + v)* J , (13) 

for j = 1, ■ ■ ■ , n, and for (t,x) £ (0, T), where fi C R n . Here the c? k are 
coupled coefficient functions which are defined explicitly via recursion. For 
each j the coefficients c£ will be defined recursively in terms of function 
c[ and their derivatives, where < l,r < k — 1. They are solutions of 
first order partial differential equations which can be solved explicitly and 
can be represented in terms of recursively defined power series under the 
assumption (|9]). We shall show that (fT3j) is valid on some domain O x (0, To]. 
Since it is desirable to have a representation which holds on the whole a 
domain 0, x (0, T] with arbitrary time T € (0,oo), in our main theorem we 
shall consider global representations of an equivalent problem, where the 
equivalence is via the time transformation r(t) : [0, 00) — > [0, 1) with 

r = (l-e~i), or t = t(r) = -/3m(l -r). (14) 

This introduces a time-dependence in the related coefficients c? k a T , even in 
the case of time- homogeneous coefficient functions x — ► b l - k {x) in @. The 
main result for parabolic systems of type ([T]) is formulated in the time- 
homogenous case, i.e. when the coefficients bj m depend only on the spatial 
coordinates. The extension to the time-dependent case is then the content 
of the subsequent corollary. 

Theorem 2.1. Given assumption ^ and some domain Q x (0, T] for any 

finite T > and any domain Q C K n there exist (3, t > such that the 
fundamental solution of 

duj ^ d 2 Uj >A . duj_ 

Q T 1 _ T 2^ dx 2 + 1 _ T 2^ °Jk d K^) 

3=1 3 j,k=l 

equivalent to |7p ( or ) via has the pointwise valid representation 

pf T (r,x,0,y) = - / J==^ex I > (- ) exp ^^.r^^^^j ' 

(16) 
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for j = 1, •■■ ,n, and for for (t(r),x) G (0,T) x 0, i.e. r G (0, 1 - e z 3 ), 
where £1 C M™. For i/ie coefficient functions ci the following holds: fork = 
we /iai>e 



cJfl T (r,a;,y) = cj(x,y) = V(?/ m - x m ) / V ft/ m (y + s(x - y))ds, (17) 
and /or all k > 1 we have 

ri l _ T 

4,/3,r( r ' x ' y) = y o K-i^A 1 ' y + s ( x ~ y)> y) sk ^ ( 18 ) 



with 



'Ezm^m( 3 



1 <9ic m L fc-1 



More explicitly, we have 

4,pA T > x >y} =4( x >y) = -Ei, m ^ b Lj(y) AxJ+h TTh\ 

= E 7 ^ 

and, given the power series representation 

4-lAr( T ' I 'S') = *|fc-l)7j(f) A ^ 7r ' 

u;e Ziaue 

iftT^ 1 .!/) = E 7 ^C(fc_i h i(y) A ^+ 

E 7 { Si Ep+a= 7 (^ + + 1 ) c ^(/3+l i ) c (A;-l-r-)(a+li) 

+ Ei(7i + 2)(7i + l)Cfc( 7+2l ) + Ep+ a=7 (E (y)x 

(ai + i)c (fc _ 1)(Q+li) } (ELopI1 Ax<5 ) » 

n 

where with 8^ := ^2 5%, &nd 

(3 



i=l 



5=0 



5=0 



(1 - T)5 S + k 



X-y-s) 



Ax d 



(19) 



(20) 



(21) 



(22) 



(23) 
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Remark 2.2. In (I2ip the notation 

<at(*> y) = E i-i )7 ^>) Aa;V ( 24 ) 

may be expected, but we dropped the /3, r indices in order to keep some 
notational simplicity in (j22l . 

Remark 2.3. Note that c? p T does not depend on r (and (3). This means 
that in (19) -j^Ck-i differs from only for k > 2. 

Corollary 2.4. Consider the same situation as in the preceding theorem, 
except that there are time dependent coefficient functions (t,x) — > b l - k {t,x) 
and in addition holds. Then a analogous statement as in the preceding 
theorem holds with recursive . 



3 Formal computation of solution of parabolic sys- 
tems of type (CO) 



First we consider the equation (TjQ) (or, equivalently, (j3])) without the time 
transformation (I14p . and with time-homogenous coefficients, i.e. where the 
coefficient functions x — ► b l - k {x) depend only on the spatial variable x. We 
consider the ansatz 



Pj(t,x,y) = * n exp 
V47rt 



■^^ + J24^y)t k . (25) 

fc=0 / 



We derive recursive relations for the coefficients cj.. In a second step, as- 
suming (J7|), and therefore global analyticity of the bj k , we derive the explicit 
solution in terms of Taylor power series of bj k . For the time derivative we 
get 



It + ^4^ + £ y^' 1 ) p &> x > y) - 

k J 



^(t,x)= ^-- + ^=f± + y 'k4(x,y)t^ \ Pj (t,x,y). (26) 
For the first and second spatial derivatives we get 



dpj _ I -Axi 
dxi ~ [ 2t 



+E^4(^)* fc W>^)> ( 27 ) 

k / 



and 



+ 



(-^ + E fe ^4(^y)* fe ) 2 V,(*^>y)- 



(28) 
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Plugging into (jl]) and ordering with respect to the terms t 2 ,t 1 etc. we get 
the following recursive relations for the where 1 < j ' < n: 

2 . Yd ^ _ V M (29 ) 

1 ■ 4t 2 ~ 2^ 4t 2 W 



- 1 -S = -E^-i(E^-E^)^). ( 3 °) 



2t ^ 2t 2t\^ dxi 

l \ l lm 

and for all k — 1 > 0. 



* 1; ^4 + ax* — ^k-l + Z~2?=lZ~2r=0\ K dxi clr dx l C k-l- 

(31) 

Note that the first order coupling of the system is essentially reflected in the 
recursive first order partial differential equations starting from (|31|) . This 
would be different if we had coupling via the second order terms and it makes 
the solution of the system much easier. Note that equation (|29p is satisfied. 
Equation (|30l) is equivalent to 

f)r> 

E A ^ = -£^>)Ax m , (32) 

l Lm 



with the solution 



4( x > y) = ~ x m) f E b im(y + s ( x - y)) ds 

m ^° I 



(33) 



and for all k > 1 we have 

-I 



cL(x,y)= I R k ^ 1 {y + s{x-y),y)s k ds (34) 



o 



with Rk-\ as in equation (|56p . Next we compute the solution explicitly 
doing the integral for cq first. We abbreviate Ax = (x — y) with components 
Axi = (x — y)i and for a multiindex a = {a\, ■ ■ ■ ,a n ) we write Ax a := 
Il" =1 Aj;" i . Furthermore, we define \a\ = Yi a i ^ 

6 L(-) = E^ m , 7 (2/)(A^, (35) 

7 
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along with some multiindex 7, then 

4(x, y) = - Em Ax m Jo Ei H m (v + sAx)ds 

= - Ei, m Ax m lo E 7 tf mi (.y)( sAx V ds 
= - Ei, m Ax m E 7 &L>)a^ Jo shlds 



E, m E 7 ^ m7 (2/)A^+ 1 ^«w+ l1 







= -E^E 7 ^ m7 (y)A^ +1 ^ 

= E 7 4yA^- 
Next we compute ci for k > 1. We have 



40^) = /o{EiE?=o 



dx; dxi 



+A Cfc _ 1 + e, &L%rK» + - v))**- 1 *- 

Assuming that cj[— 1 equals its Taylor series for every y £ R n , i.e. 

4-1^) =^ C (fc-l)7 (2/)Ax7 ' 

7 

then we may evaluate the derivatives occurring in Rk-i as follows 
<i?4-i 

-^7" = 2jrii + l)c (fc _i ){7+ i i) (y)Ax 7 , 

1 7 

and 

g^" 1 = Y^ 1 + 2 )( 7i + 1 ) C fc(7+2 I )(y)Ax 7 , 

and 

Qx~ Qx- ' = Y \ Y ^ + 1 )( ai + 1 )4( ) 8+l i ) C (fe-l-r-)(a+l i ) 
1 1 7 {l3+a=j 
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For the multiindex 7, we have 



PZ(x,y) := f (y + six-yW^da 
Jo 



(42) 



17 / n 

En 



7« 



5=0 



L - 6i)\ 



y, 



7 1 



8=0 

7 



n 



,(7-5) 



5=0 



where 5s := and s 5 = IT™ =1 s <5i = s Ss . Hence 



(43) 



c£( x >y) = 

S 7 { Si Hp+ a =y(0i + 1 )( a * + 1 ) c v( ( g+i i ) c (fc-i- r )( Q +i i ) 
+ Ei(7i + 2)(7< + l)c fc(7+2i) + E/3+a= 7 (E jtMm^y)* 

(ai + i)c( fc _i )(a+ i.)} (ELo^l Ax5 ) • 



4 Proof of convergence of the formal solution (1251) 
for some time < t < Tq 

In this section we shall show that the representation (|25D of the solution of 
holds for some time < t < Tq for some Tq > 0. In the next section 
then we will show how the time transformation t — ► r(i) is used to get a 
representation of the fundamental solution for any domain with any finite 
time horizon. We shall prove that the representation (I25|) holds for the 
equation 

dun d 2 Uj 



i^E 



1 

J= i j 



(44) 



with some (5 such that i = t/3. Essentially this step means that we have 
proved the validity of the representation for smaller time. Considering the 
solution of (j3|) at time t is equivalent to considering the solution of (|44p at 
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time r with t = t/3. If we want coefficients to be small then we choose (5 
small. Hence if 

1 / V" At 2 00 \ 

pj(r, x, y) = exp f ^ i + £ (45) 

is a representation of the fundamental solution of (|44p on some domain 
$1 x (0, T] for some r < 1 (which may be shown by proving that for each x, y £ 
f2 c? k Jx,y) converges to zero), then this is a proof that the representation 
(|25l) of the fundamental solution of (jlj) converges for t = (3t. 



4.1 Proof of convergence on bounded domains 

Since f2 C M. n is bounded, there is a ball -Br(O) around with radius R such 
that Q C Br(0). Recall that 



4( x > v) = yZiVm -x m ) f V (y + a(z - y))ds, 



(46) 



and for all k > 1 we have 

4(s,l0= / i? fe -i(y + s(a;-y),y)s fe - 1 ds (47) 
j 

with 

fl^fcz,*) = Ac^ + Er=i £to (^^4-1- 

+ Sim Hm( X )dx^ C k-l- 



(48) 



If the modulus of the coefficients bj are bounded by the generic C > for 
all j,l,m, then we have 

|cg(x, y)| <?i 2 iia (49) 
Next the time transformation 

t = (3t (50) 

transforms the equation 

du % _^d 2 ui ^ duj 

-fr-2^-d^ + l^ b 3k-for k ( 51 ) 

3=1 3 jk 



into the equation 
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where Ui(t,x) = V{(t,x), where ^ = ^rf^ = The analogous repre- 



sentation of the solution is of the form 



J* 



tfft y) = ex p ( E ^f Xi ) ( 1 + £ ) > ( 53 ) 



for j = 1, • • • , re. Plugging (|53|) into (|52|) and collecting the terms with r~ 2 , 
t -1 etc. we get (we feel free to write t instead of /3t if convenient) 

r- 2 - 2iM = ^_^. (54) 
4/3r 2 P ^4/? 2 r 2 1 ; 



*" = - £ -» E £ - £ (E ^ - E • («) 

and for all k — 1 > 

(56) 

We divide equation (|56p by /3 and get the solutions (the solution for ^ 
equals exactly that for Cq in PS1 ) 

c kB( x >y)= R{^i{y + s{x - y),y)sT*~ 1 ds. (57) 

JO 

Next we prove 

Theorem 4.1. There exists (3 > suc/i i/tai 

/or eac/i x, y G f2, 1 < j < n c? k p(x, y) j as k j oo. (58) 
Proof. First we remark that 

aw 

for some generic C, where 



;C (x,y)<cH c 7, (59) 



Cq P ■■= sup c (x,y). (60) 



Indeed, if we define 

b up := sup 6f m (x), (61) 



xBQ, l<l,m<n 
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then writing the multivariate derivative of order a with a = (a\, ■ ■ ■ ,«„), 
and \a\ := E[!=i a « we S e t \ a \ n terms bounded by C^ a){ ~ l b ulp and n 2 terms 
bounded by RC\ a % up . 

Next, a majorant of c^ k p{x,y) is obtained as follows: we consider three 

types of operators 0,' n , O k ,n , O k ,n with positive integers k, and acting on 
a single function / :QxO->Rorona families of functions (/z)i<;<fc : 
f2 x f2 — > K, namely 

Oj' n [/] (*,!/) := |A/(x,y) 

Of [/*, ...,/!] (x, y) := £ Er=i Eto f£%f ( 62 ) 



Let 
For 

we have 
where 



M k :={(a k ,--- , ai )\ aj e {1,2,3}} (63) 



c t P /3 := SU P ^pfav) ( 64 ) 
a;,j/eOj'g{l,— ,n} 



c k P p - O k' nc o( x ^)= Y O^' n c (x,y), (65) 



0^ n [f] (x,y) := O^Ofrr o • • • o 0?> n [/] (x,y). (66) 

First let (resp. 2^,3^) the multiindex a E Af^ such that for each 1 < 
m < k a m = 1 (resp. a m = 2,a m = 3). Hence 

G\ [f] (x,y) = (O 1 )" [/] (x,y) = A k [f] (x,y) (67) 
etc.. Then for all x,y £ 17 

ok k r 2k up 

|Qfc' n N < r, , (68) 



and for 6 bounded by a generic C on the domain we also have 

|o;' n N(*,y)l< * fe , ■ (69) 

The operators of quadratic type applied to co(x,y) O k ,n co decrease to zero 
as k | oo if (5 is small. We estimate (a rough estimate is sufficient here) that 
surely we have 



,2,n m , m ^ /3 fc A; 2 2 fc - 1 n^C fc A;!(c7) fc+1 

"fcf 



|0fc2 ,n [/](X)2/)| < y ^ , . (?0) 
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For large k this is essentially the largest term of all the 3 k contributions in 
the sum ()65[) for large k {k fixed). We mean the following: if we choose 



/3< 3-4n 2 C 2 (c^) 2 ' (71) 
then surely we have for k > ko (some ko > 

ok okL.'2')k-l r ,k(~iku\( r u P\k+l 

\3 k Ol [c ] (x,y)\ < 3fjk2 fc - (C ° } 4 0. (72) 

as k | oo, and this is also the estimate which holds for c& for large k. Here 
we choose j3 such that in a summand in O k ' n co(x, y) in ([65]) each occurrence 
of an operator of form Ot' n can be replaced by an operator of form O k ' n 
in order to get a majorant estimation. So in the sum in (|65p it suffices 



to concentrate on the summands consisting of concatenations of operators 
of form Ofc' n and OJ n '. For natural numbers I let us define an increasing 
sequence of numbers k\ < ki < ■ ■ ■ < ki < ki+± ■ ■ ■ , and operators 

(73) 



Ol.n /nl.w 
, , := (J, o • • • o (J , 
k l+1 k t k l+1 k t 



02,n /^i2.n ^2,11 

, , := (J, o • • • o (J , 
h+lh k l+1 k t 



Then in the summands o (I65p we have to consider the asymptotic behavior 
of values of family of operators of form 

h h O-'-oO,', o(j, . (74) 

ki+iki fe;fei-i fc3fe2 V / 

or of form 

Ofcjfc,,, o • • • o o fc ; fc2 o o k2ki (75) 

applied to co(x, y) as fc | oo. If there is only a finite occurrence of operators 
of form O k ' n in such a family ( (I74p or (|75jl ). then the asymptotic behavior 
is clearly the same as for O k co(x,y). If on the other hand there are infinite 
occurrences of operators of form O k ' n in ( (|74j) of (I75p ). then for large k 
Ofco(x,y) becomes a majorant of such a summand. Hence, the estimate 
(|72p is a majorant for large A: and proves the convergence of the series in 

(ESI). □ 



4.2 Remark on unbounded domains 

It is not possible to extend the proof in the preceding section immediately 
to unbounded domains f2 C M n . However, a similar proof with a different 
but equivalent ansatz 

lj(t,x,y) = ^exp (- £ ^ t Az? ) (l + fXfoy)**), (76) 



fc=0 
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leads to such an extension. The recursion equation for do and cq are equiv- 
alent, but the recursion equations for the dk,k > 1 are more involved. How- 
ever, it can be shown that given t, y the supremum in x of each 

1 f-Ax 2 \ • 

?r exp rr~ f ) ( 77 ) 



V4^t V 4 * 

is in some ball which can be chosen a priori. However, this is beyond the 
scope of the present paper, and we shall consider a similar situation in [7]. 

5 Extension to the time-inhomogeneous case (so- 
lution and global convergence) 

In a second step we use a certain nonlinear time transformation in order to 
lift the result to convergence for any finite time t. This requires the extension 
of the analysis to the case with time-dependent coefficients. Note that in the 
extension of the recursion of the to the time-inhomogeneous case only first 
order time derivatives occur. This is the reason for the weaker constraints for 
([5]). We start this Section with the computation of the recursive coefficients 
Cfc in the case of time- and space-dependent drift coefficients l? kl . Then we 
shall complete the proof for convergence on bounded domains for any finite 
time in the time-homogenous case, and finally in the time-inhomogeneous 
case in the following subsections. 

5.1 Formal computation of recursive coefficients in the time- 
inhomogeneous case 

We consider parabolic equations with time-dependent coefficients of the form 

jk 

We consider the ansatz 

1 / Ax 2 °° \ 

Pj (t,x,0,y) = ^=nexp -— - + ^4(t,z,y)i fe . (79) 



Compared to the time-homogenous case the time derivative contains an 
additional term. We have 



x, y) - ^ - ft + ^4*2^ + Sfclo ijtK*' x ' y) tk 



(80) 

+ J2 k kc j k {t,x,y)t k ~ 1 )p :j {t,x,0,y) 
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The spatial derivatives are essentially the same as in the time-homogenous 
case. We compute 

^ (t ' s,y)= (^ + S^(^»i/)^i(^.o J i/). ( 81 ) 

and 



+ (-^ + E k £4(t, x, y)t k ) Pj (t, x, 0, y) 



dxf \ ' > a I \ 2t 1 ^k d x 2 

(82) 

dxi 

Plugging into [Hand ordering with respect to the terms t _2 ,t _1 etc. we get 
the following recursive relations for the ci, where 1 < j < n. First, the 
highest order terms are the same as before: 

, Y\ 4 Ax, 2 v-^ Ax? , N 

The terms of order t~ l are essentially as before (we just have to add the 
t-argument in the coefficient functions b l - k ): 

Z \ i im / 

For k — 1 > we get an additional i-derivative on the right side: 

. d < _ 9c fe-i i a J i ( a J d 



l ■ KC- k -+- 2_^l tlXl g Xl — dt "I" ^Cfc_i "I- 2^1=1 2^r=0 {dx^ dx^k-l 

+ E lm HJt,x)^c{_ 1 ^ Ri^(x,y) 
Hence, 



(85) 



dc 1 

E A ^ = -E^mfo*)^, (86) 



1 (,m 



which has the solution 



CoO^y) = ~ Xm ) / ^tf mit'V + s ( x - y)) ds 

Jo , ' 



(87) 



and for all k > 1 we have 



<4( a; '2/) = / Rk-i(t,y + s(x - y),y)s k ds 
Jo 
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with Rk-i as in equation (|85p . The explicit calculation of the solution is 
know completely analogous, so it suffices to write down the results. We write 

7 

along with some multiindex 7. Then 

4&*,v) = - Ei, m E 7 b i mj (y) Axl+u TTH 

(90) 

Given that cjL-. equals its Taylor series for every y S M n , i.e. 

4-i& *)=£ i-ih(*> y ) Ax7 = E i-D>) Aa ^ ( 91 ) 
7 7,^ 

we have 

4(t> x >y) = E 7 ,« /c ( fc _i) 7 i(y) A ^ 

+ E 7 { Ei E/3+a= T (A + !)(«* + ^^OS+lijC*' ^cjfc-x-rJCo+li)^' y ) 

(92) 

+ E 4 (7 4 + 2)(7i + i)c fc(7 + 2l ) + E/3+ a=7 (E ?j[^(^y)x 
(oi + i)c( fc _i) (a+ll )} (ELo^5 Aa;5 ) . 

where the are defined exactly as before. 

5.2 Completion of convergence proof for finite time in the 
general case on bounded domains 

We apply a nonlinear time transformation. Consider for f3 > the transfor- 
mation r(i) : [0,oo) — ► [0, 1) 

T =(l-e~i), or t = t(r) = -/?ln(l - r) (93) 

with 



<9r 1 — r 
The transformation of the equation 



(94) 
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then is 

Let us call the associated coefficients of the global expansion of the fun- 
damental solution by cj, ~ T . If we can show that for each given x,y (in Q 

and then in W 1 in general) c? k g T (r,x,y) converges to zero as k | oo, then 
we have convergence for r < 1 which implies convergence of the analytic 
representation for the original equation for t £ (0, oo). First we derive the 
recursive relations for (I96p . Since 

4>i ' T (t, x, y) := 1 exp f- — (97) 



V47r(-^ln(l-r)r V 4 (-/31n(l - r)) 
is the fundamental solution of the equation 



dr 1 — r 

we consider the ansatz 



Au, (98) 



pf ' r (r, x, „) = <^(r, x, 2/) exp ]T 4,/^ v)t* ■ (99) 

\fc=o / 

Using t = — /31n(l — r) we have 

8pi( f ™ ,.\ _ [ _ n. at , Ej Ai? at ,poo j / n t; 



(100) 



and 



+ (-^ + A4j>,r(^ 2 J rf ,T fc ^ y)- 

(102) 

Plugging into (|96jl and ordering with respect to the terms r 2 , r 1 etc. leads 
to 

„- 2 . EiAxf dt _ (3 ^Ax ; 2 



4i 2 8t 



T^r E^> (^) 
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which is satisfied because the second order diffusion term in ([96]) is . For 
the terms of order r _1 we get: 



--i • _ n 9t - V A. 

2t dr 1-t l^l 2t 



(104) 



For fc — 1 > we get an additional r-derivative on the right side: 
T k-i . , i | £ V a Ti H _ 9c I-i , £ A j 

1^1=1 2^r=0 ydxi C rd Xl C k-l-rJ V 7 



We have 



Co,p, T ( T ,x,y) = yZ(ym-Xm) [ V &I |Tn (t(r), y + a(ac - y))ds, (106) 

and for all k > 1 we have 

4 AT (x,y)= f B i l fL 1 (t(r),y + s{x-y),y)s iL T }Jl - 1 d S , (107) 



where 

1 — r d ■ 

Rll^x^y) = — —— c \_ 1At + ^ T _ 1AT (T,x,y) (108) 
with R l k _i is as in equation (|105p . From this and the preceding sections it is 



clear how we get the power series representation (|22|) in theorem 2.1. above. 
We see from this representation that the proof for small t can be used, only 
that the substitution 

P > ~7~~ (109) 
1 — r 

has to be made. Since there are only first order time derivatives in the 
recursion (cf. ([85]) and A1Q5 j> . the proof of section 5.1. can be extended triv- 
ially. Hence, global convergence (for any positive t of our analytic expansion 
follows from the following 

Proposition 5.1. For each constant c > the range of the function 

{p, T )->t = -p\n(l-T) (110) 
is unbounded on the domain 

(p, T )\JL = c\. an) 
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Proof. c=j^ = -e/31n(e(l - r)) | oo as e | 0. □ 

This means that it suffices to prove that the recursion (|1U41 (|105|) con- 
verges to zero for some (which may be as small as we need). 

6 Representations of initial boundary value prob- 
lems of first and second type 

The explicit fundamental solution leads to representations of solutions for 
initial-boundary problems of parabolic systems and parabolic equations. We 
consider two examples. 

6.1 Representation of the solution for initial-boundary prob- 
lems for parabolic systems of first type 

For the Cauchy problem (|lip we have the following representation of the 
solution u: 

u(t, x) = J Rn (p(y)p(t, x; 0, y)<j>{y)dy 

(112) 

+ Jo* Jr" /( s > y)p(*> x > s ' v)4>{y)dyds 

Remark 6.1. Strictly speaking, the solution for p presented here is on bounded 
domains Q, C M. n (which is for large O a sufficient approximation for nu- 
merical purposes, but not exact). However, such exact representations on 
unbounded domains can be found using the recursion indicated in Section 
4.2. 

6.2 Representation of the solution for initial-boundary prob- 
lems for parabolic equations of second type 

In the case of the scalar problem [12] for the solution u the ansatz for 7 with 
u(t, x) = f n 4>(y)p(t, x, 0, y)dy - f Q f n f(s, y)p{t, x, 0, y)dyds 

(113) 

+ Jo Ib K*> x ; s ' Vh(s, y)dsdy 
leads to the integral equation 

±7(f, x) = J * J B I §^(t, x; s, y) + a(t, x)p(t, x; s, y) \ j(t, x)dB x ds 

(114) 

+h(t,x) 
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where 

h (t> x ) = In ^ s ' y)<t>(y)dy 
+ Jo /n x ; s > y)/( s ' 

+a(t,x) f n p(t,x;s,y)(j){y)dy ( 115 ) 

-a(t, x) Jo J n p(t, x; s, y)f(s, y)dyds 

-ip(t,x) 

Hence with our explicit solution for p we reduce the initial-boundary value 
problem of second type to the solution of a linear integral equation. 



7 Generalizations, applications, and further com- 
ments 

The preceding results can be extended to more general diffusions. We have 

Theorem 7.1. Consider equation with space- dependent diffusion coef- 
ficients x — ► a l - k {x) which satisfy 

\d» k \ < c |a| (116) 

Assume that the conditions of theorem 1 are satisfied. Then the fundamental 
solution has the representation 



P l {t,x,y) = 




(117) 



where for each i (x,y) — ► d 2 (x,y) are functionals which assign to each pair 
of points x,y the length of a geodesic with respect to the line element 

dsj = ^2g l jk dxjdx k , (118) 

jk 

with {g l j k ) the inverse of {a^ k ), and the c k are smooth functions given by 
recursive relations similar to that in theorem 1 but involving df and partial 
derivatives of d 2 . 

The proof is quite analogous except that additional existence and reg- 
ularity results for the Riemmanian metric functional d 2 are needed. These 
are given in [S]. General analytical formulas are not available for the func- 
tional d 2 but in [5] that solutions can be approximated in arbitrarily strong 
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Sobolev norms. This may be used to obtain approximations of (|117j) in ar- 
bitrarily strong Sobolev norms when combined with the results in [6]. Note, 
however, that an extension is far from obvious if the second order terms are 
coupled. 

An immediate application of theorem 7.1. is a result of Varadhan which 
we state and prove here in the case of time- homogeneous coefficients and for 
scalar equations, where the highest order coefficient function in the global 
expansion may be denoted by d 2 without an index i. 

Corollary 7.2. (time-homogeneous and scalar case) Assume that for each 
i we have A£ 2 < a % jk(x)£,i£,j < A£ 2 for x G Q C R n and some constants 
< A < A. Then for bounded Holder- continuous coefficient functions x — > 



where d 2 is the Riemannian metric induced by the line element H118]) . 

Proof. The reason for the assumption of Holder continuity is just for the 
existence of the fundamental solution, which may then be ensured by the 
parametrix method). For the assumptions of theorem 7.1 this follows directly 
from the representation (|117f) . Given x,y one may define in a bounded 
domain x,y £ Q containing the geodesic a series of coefficient functions 
( a jfc )« an d (b l j'k)n satisfying the assumptions of theorem 7.1. and such 
that a-^ix) — > a l j k (x) and (b l - k ) n — ► b % - k . Here we can assume that the 
corresponding geodesies connecting x and y are in 0, □ 

There is a deep difference between the representations considered here 
with leading terms of the form 



and direct Taylor expansions of the solution. Indeed, in [2] we saw that 
for the characteristic function (the Fourier transform of the fundamental 
solution with respect to the parameter y), where a direct Taylor approach 
seems natural, it seems that convergence results can be obtained only if co- 
efficients are of linear spatial dependence. We note that results of myself for 
scalar equations cited in [3] cannot be directly generalized to the systems 
considered here. Our results may be used to generalize the results in pQ and 
construct efficient computation schemes for related reaction diffusion equa- 
tions. A second application may be the definition of generalized Brownian 
motions (cf. [15]). This was attempted in [12] in the context of elasticity 
and the Lame equation, but not in a rigorous way. Note that Lame equation 
has coupling of second order terms, so the generalized Brownian motions as- 
sociated to (|117p would not cover these examples from elasticity (because we 



°tjk( x )> x -> b jkW 



lim — 4t In pit, x, y) = d 2 

no 



(x,y) 



(119) 




(120) 
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have no second order coupling in (I117p ). However, the functional analytic 
procedure to introduce processes as measures on path spaces using Riesz rep- 
resentation theorem leads to a new class of processes. In the special case of 
higher dimensional scalar equations expansions of the type considered here 
in a probabilistic setting have been found to be very competitive (cf. [3]). 
The results presented here are also a first step to get into deeper analysis 
of quasilinear parabolic systems, and numerical methods considered in ([BJ) 
and ( [5] ) may be extended and used together with analytical results in ( [8] ) 
and ([S]) to obtain efficient and accurate schemes for quasilinear systems. 



8 Two examples of application: multidimensional 
Burgers system and Pauli equation 

The representation of solutions for parabolic systems considered in this paper 
has a wide range of applications. One simple example is the multidimen- 
sional viscous Burgers equation with forcing. It is of the form 

— -^V 2 v+(vV)u = -VF(i,x), (121) 
at 

where x £ Q C W 1 is some domain, F is some outer force, and v de- 
scribes some velocity field. If the initial condition are of potential form, i.e. 
v(0,x) = — V$(0, x) for some potential then the velocity field remains a 
potential as time goes by and the dynamics id governed by 

d t <S> - vV 2 <£> - ^V 2 $ = F. (122) 

Applying Cole-Hopf transformation ip = 2v\nip{t,x) we get an (imaginary- 
time) Schrodinger equation of form 

d t if) - iA7 2 V> - ^"-Pty- ( 123 ) 

It is clear that our main results can be easily extended to this case with 
potential term. 

Another example of interest is Pauli' s equation which describes the 
daynamics of on electron (with spin) in the presence of a (possibly space 
and time-dependent) magnetic field: 



m dt \ifc(t,x)J 



(124) 



(ipi(t,x)\ 
W> 2 (t,z)/' 
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where L is the angular momentum B the magnetic field (possibly dependent 
on space and time), and a is the vector of Pauli's spin matrices. Furthermore, 
as usual m e and e denote the mass of the electron and the electric charge 
respectively, h is the normalized Planck constant, and 



^dia = ^ 2 sm(a) (125) 



with a = Z(x,B). Again a slight extension of our proof of Theorem 2.1 
leads to a global analytic expansion of our main result. 
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